Decoherence in a Cooper pair shuttle 
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We examine decoherence effects in the Josephson current of a Cooper pair shuttie. Dephasing 
due to gate voltage fluctuations can either suppress or enhance the critical current and also change 
its sign. The current noise spectrum displays a peak at the Josephson coupling energy and shows a 
phase dependence. 
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The Josephson effect^ consists in a dissipation-less cur- 
rent between two superconducting electrodes connected 
through a weak link^^. The origin of the effect stems 
from the macroscopic coherence of the superconducting 
condensate. Since its discovery in 1962, the research on 
devices based on the Josephson effect has been achieving 
a number of important breakthroughs both in pure^ and 
applied physics^. One of the most recent exciting de- 
velopments is probably the implementation of supercon- 
ducting nano-circuits for quantum information process- 
ing^, which requires the ability to coherently manipulate 
these devices. By now, this has been shown in several 
experiments in systems of small Josephson junctions^. 

Very recently, Gorelik et a/£>£ proposed a very appeal- 
ing setup, the Cooper pair shuttle, able to create and 
maintain phase coherence between two distant supercon- 
ductors. In its simplest realization, shown in FigfTJ the 
system is made up of a superconducting grain, externally 
forced to move periodically between two superconducting 
electrodes. Despite the fact that the grain is in contact 
with only one lead at a time, the shuttle does not only 
carry charge, as in the normal metal case^*^, but it 
also establishes phase coherence between the supercon- 
ductors. 

Aim of this work is to analyze how the presence of the 
environment affects the coherent transport in the Cooper 
pair shuttle. The intcrplayii between the periodic driv- 
ing and the environmental dephasing leads to several in- 
teresting results. By increasing the coupling to the envi- 
ronment it may result in an enhancement of the super- 
current as well as in a change of its sign (7r-j unction). In 
the last part of this Letter we propose an effective imple- 
mentation of the shuttle mechanism where the switching 
of the Josephson couplings is controlled by an external 
magnetic field. The shuttle consists of a small super- 
conducting island coupled to two macroscopic leads and 
forced to change its position periodically in time, with 
period T, from the Right (i?) to the Left (L) electrode 
and back (see Fig[T|). The grain is small enough so that 
charging effects are important, while the two leads are 
macroscopic and have definite phases 4>l,r- The moving 
island is described by the Hamiltonian 
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where Ec is the charging energy, Ej ' ' is the Joseph- 
son coupling to the left or right lead respectively, and 
n g is the gate charge. The variable h is the number of 
excess Cooper pairs in the grain and (p is its conjugate 
phase, [n, ip] = —i. The system operates in the Coulomb 
blockade regime, Ej (b) < E c (b = R,L) so that only 
the two charge states {\n = 0) , \n . — 1)} arc important. 
Their relative energy (IS.Ec) is controlled by the gate 
charge n g (t). The superconducting gap is assumed to 
be the largest energy scale in the problem, so that quasi- 
particle tunnelling can be neglected. Both Ec and Ej L ' R ^ 
are time dependent: when the grain is close to one of the 
leads, the corresponding Josephson coupling is non-zero 
(with value Ej) and the two charge states are degenerate 
(positions L and R in Fig|TJ). In the intermediate region 



(position C), Ej ' = Ej ' = 0. As in Refpi we employ a 
sudden approximation (which requires a switching time 
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Figure 1: Time dependence of the Josephson and charging 
energies in the Cooper pair shuttle. The three intervals L, C 
and R, within the period T = 2tj + 2tc, correspond to the 
situations: L) represents Ej L = Ej, n g = 1/2, 
(interaction time); C) represents Ej^ =0, n g =0, 
(free evolution time); R) Ej (L) = 0, n g = 1/2, E JR = Ej - 
(interaction time). On the left hand side, the corresponding 
position of the shuttle with respect to the leads is shown for 
each time interval. 
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in each region (see Fig.l). We further assume, as in Ref£, 
that the system is at the charge degeneracy as long as the 
island is in contact with one of the electrodes. In the in- 
termediate region (C) it is not necessary to specify the 
exact variation of n g (t), only the time integral of the en- 
ergy difference between the two charge states will enter 
in the results. 

The shuttle is coupled via the charge operator n to an 
environment described by the Caldeira-Leggett models 



Hi n t = Ai(c 



bath 
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with w\ R = tahh(Ej/T b ) (cos<p L , R , sm(j> LtR , 0) ,w c 
(0, 0, 0), 
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In Eq.([2]), Hbath is the bath Hamiltonian, with boson op- 
erators a,, aj for its i—th mode. The form of the coupling 
in Eq.@ can describe either gate voltage fluctuations^ or, 
in some limits, random switching of background charges 
in the substrate^. 

In order to analyze the transport process, we evaluate 
the time averaged supcrcurrent at steady state 



dt(i(t)) , 







and the power spectrum of the current fluctuations 

" + 00 



S(u>) = 



drS , (r)e 



where 



S(r) = -([I{t + r),I(t)\ + ) - (I(t + t)) (I{t)) . 
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In the Schrodingcr picture, the current operator is (h 
1) 



J(t) = 2eEjsm(ip-(t, L ) 6 L (t) 



(6) 



corresponding to the exchange of Cooper pairs with the 
left lead. The function 0z,(t) is defined as follows: 
0£ (t) = 1 when the grain is in the L region, and 0£ (t) = 
otherwise (the functions Q R (t) and Qc(t) are defined 
analogously). In order to evaluate Eqs.((3l HJ, we need 
to compute the reduced density matrix of the grain pit) . 
After one period the evolution of p(t) can be computed 
through a map M t defined by p(t + T) = M t p{t). The 
stationary limit is obtained by studying the fixed point of 
M.t — . Since only two charge states in the grain are rel- 
evant, the reduced density matrix can be parametrized, 
in the charge basis, as p(t) = 1/2 [1 + 3 ■ r(t)], where Gi 
(i = x, y, z) are the Pauli matrices and r,i{t) = (<Ji). The 
assumption of a stepwise varying Hamiltonian consider- 
ably simplifies the form of the map yVft, obtained as a 
composition of the time evolutions of p in the intervals 
L,C,R (see Fig. [T]). Separately for each of these intervals, 
the master equation for r{t) has the form^ 



Ht) = ]T [G k (t)r(t) 

ke{L,R,C} 



2~fj(T b )w k ] S k (t) (7) 



where the bath is taken in thermal equilibrium at tem- 
perature T&. 

Here, jj{Tb) and 7c (lb) are the temperature-dependent 
dephasing rates in the regions L, R and C, respectively, 
obtained in the Born-Markov approximation. As an ex- 
ample, for an ohmic bath with coupling to the environ- 
ment a < 1, one has^ jj(T b ) = (n/2)aEj coth(Ej/2T b ) 
and jc{T b ) = 2ixaT b . This treatment is valid provided 
that 7j.c *C T b ,Ej, and that the time interval ij(c) 
is much longer than both T^ 1 and E~^ C y In the cou- 
pling regions L and R, the only energy scale is set by the 
Josephson energy, while, during the free evolution time, 
the relevant scale is the energy difference between charge 
states. Correspondingly, all the physical quantities de- 
pend on the phases 26 = Ejtj and 2\ = J c dtAEc(t)^. 
The other important variable is the phase difference 
4> = 4>l — <Pr- The effect of damping is characterized 
by the two dimensionlcss quantities ~fjtj and "fctcr^- 
Average current- In the limiting case considered by Gore- 
lik et al&, the Josephson current does not depend on the 
dephasing rates. One expects, however, that this can- 
not be always the case. If, for example, the period T is 
much larger than the inverse dephasing rates, the shuttle 
mechanism is expected to be inefficient and the critical 
current should be strongly suppressed. In fact, we find 
a quite rich scenario, depending on the relative value of 
the various time scales and phase shifts. 

The expression of the current I(4>, 9, x, 7ji j, ictc) can 
be obtained analytically from Eqs. (|3|7| ). However, it is 
rather cumbersome and not instructive, so we prefer to 
present it in some limiting cases. A typical plot of / 
as a function of 9 and cj> is shown in Fig|2] Depending 
on the value of 9 (a similar behaviour is observed as a 
function of x), the critical current can be negative, i.e. 
the system can behave as a 7r-junction. The phase shifts 
accumulated in the time intervals L,C and R, leading to 
the current-phase relation shown in FigJ^l are affected by 
the dephasing rates in a complicated way. By changing 
"fjtj and jctc-i certain interference paths are suppressed, 
resulting in a shift of the interference pattern and ulti- 
mately in a change of the sign of the current, as shown 
in FigH 

An analysis of the critical current as a function of the 
dephasing rates reveals another interesting aspect: the 
Josephson current is a non- monotonous function of 'Jjtj, 
i.e. by increasing the damping, the Josephson current 
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Figure 2: Supercurrent (in units of e/T) as a function of 
the superconductor phase difference <j> and of the phase accu- 
mulated during the contact to one of the electrodes 9. The 
other parameters are fixed as: x — 57r/6, e ,jtj = 3/4, 



4/5. The plot is obtained for T b < Ej 




Figure 3: Average current (Tf, <C Ej) as a function of the 
dephasing rates, with cj> = — 37r/4, 9 — 7n/W, x = 5n/6. As 
a function of jjtj, the supercurrent has a not-monotonous 
behavior. Note the change of sign in the current obtained by 
varying decoherence rates in each time interval separately. 



can increase. The behavior as a function of the dephas- 
ing rates is presented in Fig|3l The presence of a maxi- 
mum Josephson current at a finite value of jjtj can be 
understood by analyzing the asymptotic behaviors in the 
strong and weak damping limits, where simple analytic 
expressions are available (in the following we do not ex- 
plicitely write the temperature dependence in jj, 7c). 
i) If the dephasing is strong, / can be expanded in powers 
of e~ ljtj and e~ lct ° and, to leading order 
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Strong dephasing is reflected in the simple (i.e. oc sin</f>) 
current-phase relationship and in the exponential sup- 
pression of the current itself. 

ii)In the opposite limit of weak damping (p/jtj <C 
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tanh 
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■y.jtj (cos0 + cos 2%) t&nd sint 



ictc 



1 + cos 4> cos 2\ 



(11) 

The current tends to zero if the coupling with the bath 
is negligible during the interaction time. In this case, in- 
deed, the time evolution in the intervals L, R is almost 
unitary, while, in the region C, pure dephasing leads to 
a suppression of the off-diagonal terms of the reduced 
density matrix p(t). As a result, in the stationary limit 
the system is described by a complete mixture with equal 
weights. The current then tends to zero in both limit- 
ing cases of large and small jjtj. Therefore one should 
expect an optimal coupling to the environment where 
the Josephson current is maximum. A regime where the 
crossover between the strong and weak damping cases 
can be described in simple terms is the limit jc ~ * 0, for 
a fixed value of 8. For example, at 8 = tt/4 the current 
reads 

2e /Ej 

I = — tanh 

T \T b 



^ 2e~~<J t J[2e- 2 "iJ t J cos</> + (1 + e~^J l ) cos 2\] sin<£ 
(1 + e- 2 T./'./)(l + e - 2 TJ*./ cos (j> cos 2\ + e- 4 TJ*J) 



(12) 



In the limit of vanishing 7,/ij, Eq. (fl"2"]) corresponds to 
the situation discussed in£. Indeed, both expressions are 
independent of the dephasing rates. The difference in the 
details of the current-phase(s) relationship are due to the 
different environment. 

In all the three cases presented here, Eqs. (|10lllll2p . 
the change of sign of the current as a function of the 
phase shifts 8 or x is evident. 

Current Noise - Cooper pair shuttling is a result of a non 
equilibrium steady state process. Therefore, to better 
characterize the transport, we analyze supercurrent fluc- 
tuations as defined in Eq.(0J. This should be contrasted 
with the standard Josephson effect^, where the supercur- 
rent is an equilibrium property of the system. We first 
consider the zero frequency noise in the two regimes of 
strong and weak dephasing. 

When the dephasing is strong, correlations on time scales 
larger than T are suppressed and the noise spectrum 
reads 



4e 2 

S(0) s trong ~ 
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2" 



-Ijt.i 



cos(20)+e- 2 ^ t V(^0,x) 
(13) 

where f(6,^,X) = cos 2 (26») - e -^ ctc cos</>cosxsin 2 (26>). 
The leading term in Eq. (fF3|) is due to the damped oscilla- 
tions in the contact regions (L,R). The phase dependent 
contribution is exponentially suppressed since it comes 
from correlations over times larger than the period. For 
weak dephasing (same limits of Eq. (|11[) ). we find 



S(0) 



4e 2 



tan 2 8 sin 2 1 
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T ■yctc 2(1 + cos</>cos2x) ' 



(14) 



which shows a much richer structure as a function of 
the phases 8 and x^- Finally, we briefly discuss the 
finite frequency spectrum in the case of strong dephasing 
(see FigU]). Superimposed to the peak at the Josephson 
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Figure 4: Current noise spectrum as a function of u> for 
T = 4tj in strong dephasing limit {^jtj = 1.2 and 'y.j/Ej — 
0.008). In the inset, we plot the spectrum in a restricted range 
of frequencies to better resolve the oscillations. 



energy, there are oscillations of frequency of the order 
of T -1 . The presence of these oscillations is related to 
the periodicity of the island motion. The modification of 
these fringes as a function of the phases is a signature of 
the coherence in the Cooper pair shuttle. 

We conclude by suggesting a possible experimental test 
of our results which does not require any mechanically 
moving part. The time dependence of the Josephson cou- 
plings and n g is regulated by a time dependent magnetic 
field and gate voltage, respectively. The setup consists 
of a superconducting nanocircuit in a uniform magnetic 
field as sketched in FigJS] By substituting the Joseph- 
son junction by SQUID loops, it is possible to control 
the Ej by tuning the applied magnetic field piercing the 
loop. The presence of three type of loops with differ- 
ent area, Al, Ah, Ac allows to achieve indipendently the 
three cases, where one of the two E/s is zero (regions 
L,R) or both of them are zero (region C), by means of a 
uniform magnetic field. If the applied field is such that 
a half- flux quantum pierces the areas Al,Ar or Ac, the 
Josephson couplings will be those of regions R,L and C, 
respectively and the Hamiltonian of the system can be ex- 
actly mapped onto that of Eq.JT]). Moreover, by choosing 
the ratios A C /A R = 0.146, and A C /A L = 0.292 the two 



Josephson coupling are equal, Ej = Ej. This imple- 
mentation has several advantages. It allows to control 
the coupling with the environment by simply varying the 
time dependence of the applied magnetic field. The time 
scale for the variation of the magnetic field should be con- 
trolled at the same level as it is done in the implementa- 
tion of Josephson nanocircuits for quantum computation 
(see Refi^ for an extensive discussion). For a quantitative 
comparison with the results described here, the magnetic 
field should vary on a time scale shorter than h/Ej, tipi- 
cally a few nanoseconds with the parameters of the first 
article in Rcf. 5 . This is possible with present day tech- 
nology^. At a qualitive level the results presented in this 
paper (7r-junction behavior, non-monotonous behavior in 
the damping) do not rely on the step-change approxima- 
tion of the Josephson couplings (which leads to Eq. (JTJ)) . 
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Figure 5: left: Sketch of the implementation of the shuttle 
process by means of a time-dependent magnetic field. Crosses 
represent Josephson junctions, right: Plot of the time vari- 
ation of the applied field (in unity of Bo = 4>o/ (2Ac),&o is 
the flux quantum) in order to realize the Cooper pair shut- 
tle. The different loop areas can be chosen in order to obtain 
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Those effects are observable even if the magnetic field 
changes on time-scales comparable or slower than Ej. 
The only strict requirement is that only one Joscphson 
coupling at the time is switched on. 
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